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We perform a new extraction for unpolarized and polarized parton distribution functions considering a flavor decom- 
positions for sea quarks and applying very recent deep inelastic scattering (DIS) and semi inclusive deep inelastic 
scattering (SIDIS) data in the fixed flavor number scheme (FFNS) framework. In the new symmetry breaking sce- 
nario the light quark and antiquark densities are extracted separately and new parametrization forms are determined for 
them. The heavy flavors contribution, including charm and bottom quarks, are also taken to be account for unpolarized 
distributions. 
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1. Introduction 

In the recent years our knowledge about the structure 
of the nucleons has improved and by the increase of 
both acceptable accuracy and the volume of data from 
deep inelastic scattering processes, new investigations 
are also in remarkable progress fljj, |2|] . 
In DIS experiments the photon transfers the electron 
vertex momentum to the proton and scatters off spin- 
pointlike quark component of it. The probability that 
the parton of flavor / carries fraction x of the struck pro- 
ton momentum is called parton distribution functions 
(PDFs) and plays a very important role to determine 
DIS cross sections. The extraction of PDFs and polar- 
ized PDFs (PPDFs) is developed to very precise QCD 
analysis in next-to-leading order (NLO) or even next-to- 
next-to-leading order (NNLO) approximation which are 
based on new model independent hypotheses lf3l— Tl3ll . 
The inability of inclusive DIS data to distinguish quarks 
from antiquarks was always the main reason of symme- 
try consideration by many theoretical groups until very 
recent years and now the growing of SIDIS for polarized 
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and Drell-Yan experiments for unpolarized data iflil - 
1911 is the motivation of considering symmetry breaking 
models in the new analysis. 

In our previous works we studied the impact of the re- 
cent very precise inclusive structure functions data from 
DIS experiments on the determination of unpolarized 
and polarized parton distributions in the standard sce- 
nario i.e. u = d = s l20l 12111 . and now we apply 
SU(2) and SU(3) symmetry breaking scenario and have 
u + d + s in determination of PDFs and PPDFs like 
what other groups have studied recently 122142511 . In the 
current analysis we apply experimental data of inclusive 
Neutral Current Deep Inelastic Scattering (NC DIS) for 
unpolarized and spin-dependant DIS and SIDIS data for 
polarized QCD fit process explained in detail in Il26ll27ll . 
The organization of the present paper is as follows: de- 
termination of unpolarized PDFs is presented in Sec. 2 
and polarized PDFs extraction is discussed in Sec. 3. Fi- 
nally in Sec. 4 we summarize and give the conclusion of 
the analysis. 

2. Determination of unpolarized parton distribu- 
tions 

The total structure function of proton F p 2 (x, Q 2 ) in MS 
factorization scheme can be written in NLO approxima- 
tion as H 

F{ NS (x,Q 2 ) + F 2 , s (x,Q 2 ) 
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+ F^ b \x,Q 2 ,ml b ), 
here the non-singlet contribution is given by 



2,NS 



'2,9 T 4?r ' 

1 + 1 + 
lg«8 + 693 



(x,Q 2 ), (2) 



and the flavor singlet contribution is 
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The contribution of heavy flavors F^ b (x, Q l ) have been 
also added in our analysis and they are taken as in 
Ref. Hit]. In the above equations a, is the strong cou- 

6(1 - z), C« = C^ s and 



pling constant, Ci 0) (z) 



the additional NLO and C!, 1 ^ are the correspond- 
ing known Wilson coefficients which can be found in 
Ref. [30]. The PDFs combinations of qt and q% and 
Q 2 ) are a ls° we ll determined in the literatures [26]. 
Here we consider symmetry breaking for u + d + s 
and a symmetry for strange sea, s — s, so our analysis 
is effected by these new assumptions. For our QCD fit 
we use the following parametrization forms of the par- 
ton distribution functions at the input initial scale Q f 
GeV 2 
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(1) and theoretical value for the n th data point of the fi 1 ex- 
periment. AN,, is known as the experimental normaliza- 
tion uncertainty and the value of N„ shows an overall 
normalization factor for the n th experiment data. In our 
global fits, we get = 1 .098 and for the total num- 
ber of used data points we put n da,a = 3279 introduced 
in Ref. JH. 

Our analysis process is accomplished using the QCD- 
PEGASUS package in the fixed-flavor number scheme 
with consideration of massless partonic flavors and 
Nf = 3 Bill . The results of fitted parton distribution 
functions, known as KKT12, and their errors at the initial 
scale are presented in Fig. Q] and regarding to the sym- 
metry breaking scenario, a comparison of our results for 
d - u and d/u as a function of x, with the results from 
other groups and experimental data is shown in Fig.|2]in 
NLO approximation. 




Figure 1: The KKT12 parton distribution functions as a function of x 
at initial scale Q^ = 2 GeV 2 in NLO approximation. 



here we take xA = x(d - u), xS = 2x(u + d + s) and 
as we mentioned above s = 5, since our used data sets 
are not sensitive to the special choice of the strange sea 
parton distributions. Due to applying some reasonable 
constraints in the parameter space of our global QCD 
fit 112611 . only 13 parameters remained free for all parton 
flavor in the final minimization. The y 2 , , al in global fit 
procedure minimization is defined as 1281 
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where n shows the number of included data points 
and £>f ata , A£>f ata , and rf eory are the value, uncertainty 
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Figure 2: Our results for d - u and d/u as a function of .v in compar- 
ison to the results from CT10 d, MSTW08 H, ABKM10 (H and 
GJR08 jH. The E866 results (tHH, scaled to Q 2 =54 GeV 2 , are 
also shown as circles. 



3. Determination of polarized parton distributions 

Generally we consider a nucleon is formed of massless 
partons that have negative and positive helicity distribu- 
tions q ± (x, Q 2 ) and the difference 



6q(x, Q 2 ) 



q + (x,Q 2 )-q-(x,Q 2 ) 



(6) 



shows how much the parton q is responsible for the orig- 
inal proton polarization and is called polarized parton 
distribution function. 

In the present analysis for PPDFs determination we 
subjoin very recent SIDIS experimental data for polar- 
ized parton densities from HERMES OH and COM- 
PASS [19] to DIS experimental data of Ref. Q since 
these additional experiments help us to apply symmetry 
breaking and recognize u, d and s separately. 
The polarized structure function g\(x, Q 2 ) is written in 
terms of a Mellin convolution of PPDFs with the rele- 
vant known Wilson coefficients AC, 
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where a s is the strong coupling constant, e q shows the 
charge of the quark flavor q and \6q, 6q, 6g} are the cor- 
responding PPDFs. For our analysis we choose follow- 
ing functional forms for polarized PDFs in the initial 
scale Ql = 4 GeV 2 

x 6u „ = N llv T] Uv x a -'' ( 1 - x) bm > ( 1 + d Ur x) , 

x6d v = N dv Tid v x a ""(l -x) h "' (l + d d x), 

x6A = N A i] A x a Hl - x) b Hl + c A Vjc), 

x 61 = N^x^il - x) b H\ + c 2 \Gc), 

x 6s = N s T] s x a °(\ - x) h '(l + d s x), 

x6g = N g T) g x a s(l-x) b s, 
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6d - 6u and 61. = 6d + 



(8) 

The nor- 
malization constants N q are determined such that the 
value of rj q become the first moments of PPDFs, i.e. 

T] q = J dx6q(x, Qq). Since the current SIDIS data are 
not sufficient yet to differ s from s , we apply 6s = 6s 
throughout, also we have to make some constraints on 
the parameter space to control the x dependance of 
PPDFs 12711 like what we do for unpolarized PDFs. 
The value of parameters tj Uv and rj^ shows the first mo- 
ments of polarized valence quark distributions 6u v and 
6d v which can be linked to F and D determined in neu- 
tron and hyperon ff-decavs 13511 by assuming S U(2) and 
SU(3) flavor symmetries 112411 . These quantities result 
into t] Uv = +0.928 ± 0.014 and i] dr = -0.342 + 0.018 
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Figure 3: The quark polarized distributions as function of x comparing 
with DSSV09 1 24] and LSS10 lH models in NLO approximation. 



as shown in Ref. M20H . Since in the present analysis we 
are not interested to force S U(2) and S U{3) flavor sym- 
metry, we should relax the symmetry relations in T] Uvzdv 
measurements by introducing two flexible parameters, 
e su(2 and e like what DSSV09 112411 has proposed 

AL u -AL d = (F + D)[l+e lvm \, (9) 

AE M + AL d - 2AE S = (3F - D) [l + e, m ] . (10) 

In above equations e su , 23 . determine the deviation value 
from S U(2) and S U{3) symmetries and are also consid- 
ered in the QCD global fit as free parameters. 
Our polarized analysis is done using the QCD- 
PEGASUS package in the fixed-flavor number scheme 
with consideration of massless partonic flavors and 
Nf = 3 same as unpolarized procedure 113 ill . Finally 



our minimization for ;^ lobal is performed with 15 un- 
known parameters from PPDFs parametrization forms 
and we obtain = 0.829 which shows an acceptable 
fit to the number of 491 experimental data. Fig. [3] shows 
the comparison of extracted PPDFs with other models 
and the symmetry breaking effect on 6u and 6d differ- 
ence, comparing with the results from other models and 
experimental data, is presented in Fig. [4] 



4. Summary and conclusions 

In the present paper we present two NLO QCD analy- 
sis of the unpolarized and polarized data from DIS and 
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Figure 4: The quark polarized distributions for the difference of x(Su- 
Sd) at Q 2 = 2.5 , 3 GeV 2 comparing to results from other models and 
experimental data 1 18, 2II2H1, 



SIDIS experiments. While the analysis we always have 
S U(2) and S U{3) symmetry breaking i.e. 11 + d + s, 
but we consider s = s since the current available ex- 
perimental data are not yet enough to recognize them. 
The effect of symmetry breaking in determining PDFs 
and PPDFs is shown and also we find out that the gluon 
helicity is still not well known 12711 . Having extracted 
PDFs and PPDFs, we can determine nucleon unpolar- 
ized and polarized structure functions Fo and g\ too. In 
general our results are in good accord with other mod- 
els determinations and this proves the progress of the 
way toward a precise description of the unpolarized and 
polarized parton component of the nucleon. 



Acknowledgements 

A.N.K thanks CERN TH-PH division for the hospital- 
ity where some parts of this work was accomplished. 
The current project was supported financially by Sem- 
nan University and the School of Particles and Accel- 
erators, Institute for Research in Fundamental Sciences 
(IPM). 



References 

[1] A. De Roeck and R. S. Thorne, Prog. Part. Nucl. Phys. 66 (201 1) 
727. 



[2 
[3 

[4 

[5 

[& 

[7 
f8 
[9 
[10 

[11 

[12 
[13 
[14 
[15 

[16 
[17 
[18 
[19 
[20 
[21 
[22 
[23 
[24 
[25 
[26 
[27 

[28 

[29 
[30 
[31 
[32 

[33 

[34 
[35 



D. Boer etal, arXiv:l 108.1713 [nucl-th]. 

W. K. Tung, H. L. Lai, A. Belyaev, J. Pumplin, D. Stump and 

C. P. Yuan, JHEP 0702 (2007) 053. 

H. L. Lai, M. Guzzi, J. Huston, Z. Li, P. M. Nadolsky, J. Pumplin 
and C. P. Yuan, Phys. Rev. D 82 (2010) 074024. 
A. D. Martin, W. J. Stirling, R. S. Thorne and G. Watt, Eur. Phys. 
J. C 64 (2009) 653. 

A. D. Martin, W. J. Stirling, R. S. Thorne and G. Watt, Eur. Phys. 
J. C 70 (2010)51. 

G. Altarelli, R. D. Ball, S. Forte and G. Ridolfi, Acta Phys. 
Polon. B 29 (1998) 1145. 

R. D. Ball, G. R idolfi, G. Altarelli and S. Forte, 
arXiv:hep-ph/9707276 

C. Bourrely, F. Buccella, O. Pisanti, P. Santorelli and J. Softer, 
Prog. Theor. Phys. 99 (1998) 1017. 

H. Khanpour, A. N. Khorramian and S. A. Tehrani, PoS EPS 
-HEP2011, (2011)443. 

H. Khanpour, A. N. Khorramian and S. Atashbar Tehrani, Few 
Body Syst. 52, (2012) 279. 

S. Atashbar Tehrani and A. N. Khorramian, JHEP 0707, (2007) 
048. 

A. N. Khorramian and S. Atashbar Tehrani, JHEP 0703, (2007) 
051. 

F. Arbabifar, A. N. Khorramian, S. Taheri Monfared and 
S. Atashbar Tehrani, Int. J. Mod. Phys. A 26, (201 1) 625. 
A. N. Khorramian, S. Atashbar Tehrani, F. Olness, S. Taheri 
Monfared and F. Arbabifar, Nucl. Phys. Proc. Suppl. 207-208, 
(2010) 65. 

R. S. Towell et al. [FNAL E866/NuSea Collaboration], Phys. 
Rev. D 64 (2001) 052002. 

J. C. Webb et al. (NuSea Collaboration), arXiv: 
hep-ex/0302019 

A. Airapetian et al. [HERMES Collaboration], Phys. Rev. D 71 
(2005) 012003. 

M. G. Alekseev etal. [COMPASS Collaboration], Phys. Lett. B 
693 (2010) 227. 

A. N. Khorramian, S. Atashbar Tehrani, S. Taheri Monfared, 
F. Arbabifar and F. I. Olness, Phys. Rev. D 83 (2011) 054017. 
H. Khanpour and A. N. Khorramian, Acta Phys. Polon. B 41 
(2010) 2929. 

D. J. Broadhurst, A. L. Kataev and C. J. Maxwell, Phys. Lett. B 
590, 76 (2004) |hep-ph/0403037]. 

A. L. Kataev, Phys. Part. Nucl. 36, S168 (2005) 
|hep-ph/0412369|. 

D. de Florian, R. Sassot, M. Stratmann and W. Vogelsang, Phys. 
Rev. D 80 (2009) 034030. 

E. Leader, A. V. Sidorov and D. B. Stamenov, Phys. Rev. D 82 
(2010) 114018. 

H. Khanpour, A. N. Khorramian and S. A. Tehrani, 
larXiv:1205.5T94 [hep-ph]. 

F. Arbabifar, A.N. Khorramian and S. Atashbar Tehrani, Pub- 
lished in the proceeding of: XIV Advanced Research Workshop 
on High Energy Spin Physics, DSPIN-1 1, (2012) 41. 

A. D. Martin, W. J. Stirling, R. S. Thorne and G. Watt, Eur. Phys. 
J. C 63 (2009) 189. 

M. Gluck, C. Pisano and E. Reya, Eur. Phys. J. C 50 (2007) 29. 
W. L. van Neerven and A. Vogt, Nucl. Phys. B 588 (2000) 345. 

A. Vogt, Comput. Phys. Commun. 170 (2005) 65. 

S. Alekhin, J. Blumlein, S. Klein and S. Moch, Phys. Rev. D 81 
(2010) 14032. 

M. Gluck, P. Jimenez-Delgado and E. Reya, Eur. Phys. J. C 53 
(2008) 355 . 

B. Lampe and E. Reya, Phys. Rept. 332 (2000) 1 . 

C. Amsler et al. (Particle Data Group), Phys. Lett. B 667 (2008) 
1. 



4 



